Proximity of Iron Pnictide Superconductors to a Quantum Tricritical Point 
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We determine the nature of the magnetic quantum critical point in the doped LaFeAsO using 
a set of constrained density functional calculations that provide ab initio coefficients for a Landau 
order parameter analysis. The system turns out to be remarkably close to a quantum tricritical 
point, where the nature of the phase transition changes from first to second order. We compare 
with the effective field theory and discuss the experimental consequences. 

PACS numbers: 71.10.Hf, 74.20.Pq, 74.70.Xa, 74.40.Kb 
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Layered FeAs materials have been extensively studied 
since the discovery of superconductivity with a transi- 
tion temperature up to 28 K in LaFeAsOi_a;F x and ex- 
ceeding 50K in related compounds. Neutron scattering 
experiments [l[ have shown that in the FeAs planes of 
the layered parent compound long-range magnetic stripe 
order develops [Fig. HJa)] which was predicted by den- 
sity functional theory [2j. Upon doping the compound 
magnetic order is suppressed and disappears at zero tem- 
perature at a quantum critical point (QCP). The phase 
diagram is strikingly similar to that of heavy fermions 
where the QCP is inside or on the edge of the super- 
conducting dome [3j. The role of a QCP in promoting 
superconductivity has been stressed in heavy fermions [4] 
and cupratesQ which renders the nature of the QCP of 
fundamental importance. 

Empirically the magnetic thermal transitions show ba- 
sically two types of behavior In some compounds, 
particularly of the 122 family, such as SrFe2As2, the tran- 
sition is first order, while in compounds of the 1111 family 
such as LaFeAsO it appears second order like. This sug- 
gests that the iron pnictides are close to a tricritical point, 
i.e. a point in the T-x plane, with x a non-thermal pa- 
rameter, where the nature of the transition changes from 
first to second order. It is in principle conceivable that 
by changing two non-thermal parameters the tricritical 
point is driven to zero temperature producing a quantum 
tricritical point (QTCP) where a range of unconventional 
quantum critical phenomena is expected to occur. Acci- 
dental proximity to such a QTCP can also dominate the 
finite temperature crossovers and is believed to occur in 
several compounds [§]. 

Whether such a scenario is viable for the doped iron 
pnictides can be determined from a Landau order pa- 
rameter analysis if the coefficients that appear in the 
expansion are known. Identifying all possible and rel- 
evant magnetic phases close to the potential tricriti- 
cal point, we have determined all the coefficients of the 
Landau theory in LaFeAsO as a function of doping from 
first principles, computing total energies in a constrained 
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FIG. 1: (Color online) Competing magnetically ordered states 
in LaOFeAs. 



density functional approach [1 01) within the local-density 
approximation (LDA) [i"i"l - ll3lj . We find that LaFeAsO 
is surprisingly close to a QTCP which will strongly af- 
fect superconducting and normal state properties. The 
resulting effective field theory [l4| reveals that at this 
critical point an Ising and a continuous order parameter 
vanish concomitantly. Upon doping a structural transi- 
tion must therefore be very close to the magnetic one, 
as indeed is found in the experimental phase diagram of 
LaFeAsO. 

The Landau theory is based on the assumption that the 
dominant magnetic instability of FeAs planes is at mo- 
mentum (w, 0) and (0, tt) where we use a notation with a 
single Fe per unit cell and take the Fe-Fe distance a = 1. 
Besides the well known magnetic stripe phase one needs 
to include here also the phase with orthogonal magnetic 
moments, termed orthomagnetic (OM), and the phase 
with spin and charge order (SCO) (c.f. Fig. [T]) [9(. In 
accord with the Landau analysis, these states appear as 
low-lying energy phases in microscopic computations (l5| . 
ab initio [l6| and density matrix renormalization stud- 
ies [rzj]. 

Within the Landau expansion the order parameter is 
given by the six Cartesian components of the Fourier 
transform of the magnetization at the relevant wave vec- 
tors. We need to keep the expansion up to the sixth 
power of the order parameter in order to account for sec- 
ond and weakly first order phase transitions close to a 
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tricritical point. Using symmetries it is found that the 
Landau expansion has 7 independent coefficients which 
we need to determine [§]. The energies of the three or- 
dered phases shown in Fig. [T] as a function of their total 
magnetization, Mt, completely determine the 7 Landau 
coefficients. The energies are given by 
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For the SCO state the factor 1/2 accounts for the fact 
that only one half of the Fe sites are magnetized. 

For the LDA computations we used the Vienna ab- 
initio simulation package (VASP) [lq] ■ The Kohn-Sham 
equations in the self-consistent calculations have been 
solved using the projector augmented wave method (l9j 
with the valence pseudo-wave-functions expanded in a 
plane wave basis set with a cut-off energy of 500 eV. All 
the integrations in the Brillouin zone are performed ini- 
tially with a Gaussian smea ring method and then checked 
with a tetrahedron scheme [20( using a sampling grid of 
10 x 10 x 6 k-points. 

For iron pnictides the method of choice is LDA be- 
cause it provides magnetic properties closer to experi- 
ment than, for example, generalized gradient approxima- 
tion [21] . We use experimental lattice and internal pa- 
rameters fixed at zero doping [l[ with symmetry group 
P4/nmm, constructing our unit cell with 4 Fe sites to be 
able to allocate the mentioned magnetic structures. We 
also performed calculations changing the z coordinate of 
the As as explained below. Electron doping has been 
introduced in our computations by the virtual crystal 
approximation [22$ . 

We fixed both the modulus and direction of the mag- 
netization to the patterns dictated by the Landau theory 
by implementing in VASP the possibility to perform LDA 
contrained calculations fiojj and computing the total en- 
ergy as a function of My. The local magnetic moments, 
needed to implement the constraint, where found by in- 
tegration of the magnetization density in atomic Wigner- 
Seitz spheres centered at the Fe sites. The energy versus 
Mt curves were fitted with expressions Eq. ^ to deter- 
mine the Landau coefficients. The minimum of the Lan- 
dau energy determines the equilibrium magnetization. 

The MS phase breaks C4 symmetry thus we expect 
that the lattice will distort as indeed observed experimen- 
tally. However we are interested in the behavior close to 
the QCP between the magnetic and non-magnetic state, 
where the orthorhombicity becomes negligible [23[ thus 
for simplicity we neglect this effect. 

Fig. [5] shows the total energy per Fe for the MS and 
OM state as a function of the local Fe magnetic moment 
Mt and different doping concentrations. The zero of the 




x= 


(b) / , 


0.25 ■ 




0.281 • 


'* ll 


0.297 » 




0.375 ' 

















* 

(c) fj 




■*/ 


"""H3^---I;.x.-j 




/= , 




■ X * — • 




0.125 • \ 




0.219 * 




0.25 - 





0.5 1 1.5 

M T (|1 B ) 




0.5 1 1.5 

M T (H B ) 



FIG. 2: (Color online) Energy per Fe for the MS (a),(b) and 
the OM phase (c),(d) as a function of the Fe magnetic mo- 
ment for different doping concentration x. The points are the 
results of LDA calculations whereas the continuous lines are 
fits using Eq. JTJ. 



energy is taken at the Mt = state. The points are 
the LDA data and the lines the Landau fits. In general 
we find that the latter provides an excellent fit to the 
data, even when Mt is not small. For the MS in the un- 
doped case (x = 0), the energy curve shows a clear deep 
minimum for Mt ~ 1-5/iB [Fig. |5Ja)]. By increasing the 
F-doping, the minimum moves to slightly lower moment 
values and eventually disappears for x > 0.3 [Fig. |2£b)]. 
The energy is extremely flat close to the critical point 
with the large moment state almost degenerate with the 
zero magnetic state and without the appearance of a no- 
ticeable barrier. This behavior already indicates proxim- 
ity to a QTCP and an anomalously "soft" magnetism at 
the transition. For the OM state Fig. (2fc),(d) the be- 
havior has a more pronounced first order character. The 
metastable minimum persists up to large dopings with a 
sizeable barrier separating the large magnetization state 
from the low magnetization state. Finally for the SCO 
state (not shown) the evolution is that of a typical second 
order phase transition. 

From the statistical mechanics point of view the LDA 
approximation is a mean-field theory Thus in the second 
order region of the phase diagram, at the LDA level, one 
finds classical critical exponents. The order parameter as 
a function of a non-thermal parameter should vanish as 
\x — x a \P with Plda = 1/2 far from the tricritical point in 
the second order region and (3lda = 1/4 at the tricritical 
point. In the critical region the energy landscape is very 
flat which is a big obstacle for the converge of conven- 
tional LDA computations. We avoid such complications 
by constrained computations that allow to extract the 
intrinsic LDA behavior. 

The susceptibility a -1 of the nonmagnetic phase at 
momentum (jr, 0) or (0, tt) should be independent of the 
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FIG. 3: (Color online) (a), (b),(c): Landau parameters for the 
three possible magnetic phases as function of the doping con- 
centration x. Dots are the numerical results while the lines 
are fits as explained in the text. The left (right) arrow indi- 
cate the point where Bi — B2 — (a = 0). Panel (d) shows 
the LDA phase diagram for the MS in the z(As)-doping plane. 



phase under consideration, thus a single parameter a ap- 
pears in Eqs. (p}. As a consistency check we allowed for 
different values of a in the fits and found that indeed a 
converges to practically the same values as a function of 
doping except for x ~ and 0.5 where higher order terms 
in the expansion become important [Fig.^a)]. The van- 
ishing of a at a critical doping x a determines the limit of 
stability of the paramagnet coming from large x. For the 
experimental lattice constants we find a ~ 0.26(x — x a ) 
eV//x| with x a = 0.27. 

Fig.[3Jb),(c) show the behavior of the quartic and sixth 
order Landau parameters as a function of the doping 
concentration x. The behavior is rather smooth and 
can be simply captured assuming a quadratic polyno- 
mial expansion except for the OM where a cubic term 
becomes important far from the critical point. For the 
MS and the OM phases we find that the quartic coeffi- 
cient B of the Landau expansion becomes negative above 
xb = 0.24 < x a - Therefore the transitions from the non- 
magnetic state to the MS and OM phase are first order 
while for the SCO state B3 > and the transition is a 
conventional second order one. One can judge the rele- 
vance of the QTCP in the thermal crossovers by comput- 
ing the height of the energy barrier at the point in which 
the magnetic and non-magnetic solution become degen- 
erate. For the MS the barrier from the fits is nominally 
~ 2K per Fe atom which is much below the limit of accu- 
racy of the computation, thus for all practical propose the 
MS-non-magnetic transition occurs at a QTCP in LDA. 
The low barrier reflects an almost vanishing metastability 
range around the transition and the physics will be dom- 
inated by QTCP behavior. The OM state has a larger 
range of metastability (c.f. Fig. @| but still with a negli- 
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FIG. 4: (Color online) Energy per Fe (a) and Mt (b) in 
the unconstrained LDA (dots) and the Landau theory (lines). 
The meaning of arrows is the same as in Fig. [3] 



gible barrier at the transition point ( ~ \QK). 

Fig.HJa) shows the energy of the different phases. The 
end of the line indicates the spinodal point, i.e. the point 
at which the state disappears as a saddle point solution of 
the Landau equations. The MS stripe is the most stable 
phase except close to the transition to the paramagnet 
where the OM phase becomes stable in a small doping 
interval. The OM preserves C4 symmetry and has the 
same structure factor as an incoherent superposition of 
(it, 0) and (0, ir) twins of the MS state, so it is difficult to 
distinguish with magnetic neutron scattering alone. It is 
interesting that for 122 compounds a state with magnetic 
order but without detectable orthorhombicity has been 
reported 24J which can be taken as a signature of the OM 
state. However this is the region where the orthorhom- 
bicity is expected to be smaller and further experimental 
work is required to clarify this point. 

Fig. 0|b) shows Mt vs. doping. Despite the abrupt 
suppression of the order parameter from My ~ 1 Ms to 
Mt = for the MS, the transition is weakly first order. 
The OM magnetization shows a sharper first order be- 
havior although, as discussed above, the barrier is very 
small. The SCO behaves as a typical second order tran- 
sition. The behavior of the order parameter of the MS 
is consistent with the sudden drop of the magnetization 
as a function of doping observed in this compound [2o| . 
Such behavior reinforces our conclusion that the system 
is close to a QTCP although our critical doping is overes- 
timated as discussed below. The fact that the OM state 
lies only less than lOmeV above the paramagnetic state 
for a large range of doping suggests that fluctuations to 
this state may be the most relevant ones in the supercon- 
ducting region. 

The LDA critical doping ~ 0.26 is larger than the ex- 
perimental one 0.05 ~ 0.06 in this compound. Since 
LDA neglects fluctuations the Landau parameters de- 
rived should be considered to be "bare" parameters. One 
well known effect of fluctuations is to reduce the stability 
of the ordered phases 26|, 27 1 shifting the critical doping 
in the correct direction. Also because of the small en- 
ergies involved, the actual transition will be sensitive to 
details as the choice of the functional or changes due to 
the relaxation of the structure. In particular it is known 
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that the magnetism is very sensitive to z(As) the z co- 
ordinate of the As in the unit cell ljl 2l|, [28| . Relaxing 
z(As) results in a decrease of around 0.1c (c = 0.87A). 
We made also computations with selected values of z and 
constructed the zero temperature phase diagram shown 
in Fig. (3^d) which shows the location of the QTCP in 
the z(As)-doping plane. We see that a decrease of z(As) 
indeed shifts the transition to the right range of doping 
moving from the weakly first order region (right of the 
QTCP) to the second order region. We also show the line 
Si = 0. We see that even if one takes a relaxed z and 
a smaller critical doping still the transition is dominated 
by the proximity to this line. Indeed we find that the 
B\ coefficient depends weakly on z{As)/c so for doping 
x = 0.05 and z/c = 0.6307 we find B x ~ 10 meV/[i%. 
Taking a characteristic moment of M* ~ 0.5/is this im- 
plies a negligible energy scale Bi(M*) i ~ 7K which again 
points to soft magnetism due to proximity to the QTCP. 

Ref. [TH has analyzed the field theory describing the 
quantum phase transition in iron pnictidcs in terms of 
fields 4>\ and 4>2 representing Neel order parameters in 
different sublattices and related to our order parame- 
ter by (pi.2 = Mi ± M2. The second and fourth or- 
der terms are precisely equivalent to ours: $3 i=1 u|</>j| 4 + 
7i|(/>i| 2 |02| 2 — ctx(4>i-4>2) 2 ■ In terms of our coefficients 
u = £>3 / 2, 71 = £>2 — B% and ax — B2 — B\. Increasing 
the doping, for the experimental atomic positions, ax 
changes from positive to negative at ib [Fig. G2h)] im- 
mediately before the transition. This is consistent with 
the appearance of the OM state in a small doping interval 
as shown above. A small ax implies a structural transi- 
tion very close to the magnetic transition as a function 
of doping 1J| which is consistent with the experimental 
phase diagram of doped LaOFeAs [25| . 

To conclude, we have shown that magnetism in an iron- 
pnictide superconductor is surprisingly close to a QTCP 
and we have determined ab initio the coefficients of a 
Landau expansion around it. The energy landscape is 
anomalously flat close to the zero temperature magnetic- 
non-magnetic transition giving rise to a very soft behav- 
ior of the order parameter in the sense that can experi- 
ence large changes as a result of weak perturbations. We 
believe frustration plays and important role in this result 
as it tends to turn Stoner like transitions into weakly first 
order ones While we have focussed on a 1111 com- 
pound one expects similar behavior to occur in 122 com- 
pounds. In particular it has been shown that magnetism 
can be quenched and superconductivity appears in the 
undoped compound by applying pressure |29|. Thus ex- 
periments exactly at the QTCP should be possible by a 
combination of pressure and doping to elucidate the effect 
of this unusual critical behavior on superconducting and 
normal state properties. The degeneracy of the ground 
state close to a QCP is believed to boost superconductiv- 
ity as a way to remove the residual entropy [il [30j| . This 



effect should be enhanced close to a QTCP, where the 
degeneracy is even larger, and may play an important 
role in determining the high critical temperature of iron 
pnictides. 
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